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ABSTRACT 

K-band spectra of young stellar candidates in four southern hemisphere clusters have 
been obtained with the near-infrared spectrograph GNIRS in Gemini South. The clus- 
ters are associated with IRAS sources that have colours characteristic of ultracompact 
HII regions. Spectral types were obtained by comparison of the observed spectra with 
those of a NIR library; the results include the spectral classification of nine massive 
stars and seven objects confirmed as background late-type stars. Two of the studied 
sources have iGband spectra compatible with those characteristic of very hot stars, 
as inferred from the presence of ClV, Nlli, and Nv emission lines at 2.078 /im , 2.116 
/im, and 2.100 /im. respectively. One of them, 116177 JRS1, has a if -band spectrum 
similar to that of Cyg OB2 7, an 03If* supergiant star. The nebular K-band spectrum 
of the associated UC Hn region shows the s-process [Knil] and [Seiv] high excitation 
emission lines, previously identified only in planetary nebula. One young stellar ob- 
ject (YSO) was found in each cluster, associated with either the main IRAS source 
or a nearby resolved MSX component, confirming the results obtained from previous 
NIR photometric surveys. The distances to the stars were derived from their spectral 
types and previously determined JHK magnitudes; they agree well with the values 
obtained from the kinematic method, except in the case of IRAS15408-5356, for which 
the spectroscopic distance is about a factor two smaller than the kinematic value. 

Key words: stars: early type - ISM: Interstellar medium: compact HII regions - 
near-infrared: young massive stars 



1 INTRODUCTION 

Massive star forming regions are commonly found embedded 
in high-density molecular clouds. They can be traced, for ex- 
ample, by the presence of CO (Churchwell et al. 1992, May 
et al. 1993, Hofner et al. 2000), CS (Churchwell, Walmsley & 
Wood 1992, Bronfman, Nyman & May 1996), and NH 3 (Ce- 
saroni, Walmsley & Churchwell 1992, Cesaroni et al. 1994), 
transitions in the radio spectra of the clouds, and are often 
associated with methanol and water masers. The presence of 
young massive stars still embedded in their parental molec- 
ular cloud can also be inferred from radio continuum and 
line surveys of IRAS sources selected according to specific 
colour criteria (Wood & Churchwell 1989; Caswell & Haynes 
1987). 
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Due to the large column density of gas and dust, young 
clusters are often affected by high visual extinction, some- 
times over 30 magnitudes, which make them undetectable 
at optical wavelengths. Dutra et al. (2003) identified many 
previously unknown stellar cluster candidates in the Near- 
Infrared (where the extinction is about a tenth of that in 
the visual window), using the Two Micron All Sky Sur- 
vey (2MASS). The 2MASS images have been used not only 
to discover new young clusters, but also to obtain physical 
properties of their members (Borissova et al. 2003, Leistra 
et al. 2005), although the limited spatial resolution (~ 2") 
restricts its usage to non-crowded regions. Several authors 
have taken advantage of the higher spatial resolution pro- 
vided by the new generation of near-infrared (NIR) array 
detectors to study clusters with high stellar surface density 
(Horner, Lada & Lada 1997, Gomes & Kenyon 2001, Han- 
son, Luhman & Rieke 2002, Massi, Lorenzetti & Giannini 
2003, Balog et al. 2004, Kumar, Kamath & Davis 2004, Lada 
& Muench 2004, Whitney et al. 2004, Figueredo et al. 2002, 
2005, Arias, Barba & Morrell 2007, Barba & Arias 2007). 





Figure 1. NIR finding charts (taken from previous NIR imaging surveys) of the clusters associated to IRAS09 149-4743 (a), IRAS15408- 
5356 (b), IRAS16132-5039 (c) and IRAS16177-5018 (d) sources. Each image is about 2' X 2.5' in size. North is to the top, cast to the 
left. In each cluster, we indicate the slit position (dotted lines) and the sources for which we obtained K band GNIRS spectra. Based on 
the GNIRS acquisition images, we found that IRAS09149-4743IRS1 and IRAS15408-5356IRS3 are in fact double, as shown in the insets 
in panels (a) and (b), respectively, at plate scale of 0.15" per pixel. 



Roman-Lopes, Abraham & Lepine (2003); Roman- 
Lopes & Abraham (2004a,b), Roman-Lopes & Abraham 
(2006a, b), Ortiz, Roman-Lopes & Abraham (2007), Roman- 
Lopes (2007), used the NIR camera CamlV, attached to the 
telescopes of the Pico dos Dias Observatory, to study clusters 
associated with high density molecular clouds, and IRAS 
sources with colours of ultracompact HII regions. In these 
studies, JHK colour-colour (C-C) and colour-magnitude (C- 
M) diagrams were used to select the cluster member candi- 
dates. Their spectral types were estimated by dereddening 
magnitudes and colours in the C-M diagram up to the point 
where they intercept the ZAMS line. This method depends 
on an accurate knowledge of at least three parameters: (1) 
the distance to the cluster; (2) the extinction law in that 
direction, and (3) the excess emission due to the possible 



presence of circumstellar material around each individual 
star, which affects mainly the Jf-band magnitudes. 

Most studies of embedded clusters assume kinematic or 
statistical distances and make use of a standard interstellar 
extinction law (Rieke & Lebofsky 1985), even though devi- 
ations have been reported in dense molecular clouds (Tapia 
1981, Indebetouw et al. 2005, Nishyama et al. 2006). In ad- 
dition to that, it is well established that massive young stel- 
lar objects (YSOs) have large infrared excess due to the 
presence of warm circumstellar dust (Grasdalen et al. 1975, 
Lada & Adams 1992), which will be reflected in an incor- 
rect estimation of their spectral types. To circumvent part of 
these problems, spectral types of some of the cluster mem- 
bers can be determined from K-h&nd spectroscopy and used 
to determine the distances to the clusters assuming different 
extinction laws. 
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In this work, we present K-band spectra of a sample of 
southern massive star candidates obtained with the Gemini 
Near Infrared Spectrograph (GNIRS). All targets were se- 
lected from previous NIR studies (Roman-Lopes et al. 2003, 
Roman-Lopes & Abraham 2004a,b, Ortiz et al. 2007). With 
the present study we intend to obtain spectral classification 
of the most massive candidates in these clusters, and from 
them the cluster distances using the spectroscopic parallax 
technique. These distances are compared to those derived 
from kinematical methods, and can be eventually used to 
improve the rotation model of our Galaxy. 

The paper is organized as follows: Section 2 presents a 
summary of the four studied regions; Section 3 reports the 
observations, data reduction and results. In Section 4 we 
discuss the results obtained for each cluster and in Section 
5 we present our conclusions. 



2 THE STELLAR CLUSTERS 

The sources chosen for this work belong to young stellar 
clusters in the southern hemisphere, which are associated 
with the IRAS09149-4743, IRAS15408-5356, IRAS16132- 
5039 and IRAS16177-5018 sources, previously studied using 
infrared imaging techniques. 

The first cluster, associated with IRAS09149-4743, be- 



longs to the Vela Molecular Ridge (VMR), and is proba- 
bly related to the optical HII region RCW 41. Ortiz et al. 
(2007) obtained JHK photometry of this cluster at 1.3" spa- 
tial resolution. The authors suggested two stars as the more 
likely candidates to ionise the nebula: IRS1, located at the 
centre of the IRAS error ellipse, and IRS2, a member of a 
small "sub-cluster" containing 6 stars, situated 1.1' south- 
cast of the IRAS position and associated to the MSX6C- 
G270. 2795+00. 8353 source. They also found another bright 
source, IRS3, the reddest object in the cluster. The K-band 
image of the region obtained by Ortiz et al. (2007) is pre- 
sented in Figure la, where these three program stars are 
identified, together with other stars that fell into the spec- 
trograph slit. 

The second cluster is associated with the IRAS source 
15408-5356 and the HII region RCW 95; it is seen against 
the Sagittarius-Carina and Scutum-Crux spiral arms. This 
cluster was studied in detail by Roman-Lopes & Abraham 
(2004b), who obtained JHK photometry of the sources in 
the region. A nebula is clearly visible in their images, with 
a clump of embedded stars that include two of the probable 
ionising sources of RCW 95, IRS1 and IRS3, as shown in 
Figure lb. 

The third stellar cluster, is located towards the IRAS 
source 16132-5039 and is associated with the HII region 
RCW 106. The infrared image obtained by Roman-Lopes & 
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Figure 3. Plots of the intensity of the Br7 and Hel line emission as function of the spatial coordinate. The emission profile is shown in 
(a), the continuum emission profile in (b) and the resulting continuum subtracted Br7 line is shown in (c). The analogous for the Hel line 
are shown in (d), (e) and (f), respectively. At the bottom we also present the contour diagram made from the GNIRS acquisition image. 
There we indicate the IRS1 and IRS7 sources, and the position of the slit together with the scale that represents the spatial coordinate 
(measured in pixels). From the analyses of this diagram we were able to estimate the amount of nebular contamination and to correct 
the stellar spectrum from it. 



Abraham (2004a) shown in Figure lc, reveals a spheroidal 
nebula containing a bright star at its centre, IRS1, and a 
smaller concentration of stars to the south-west of the IRAS 
source, which coincides with the mid-infrared source MSX5C 
G332. 5302-00. 1171; the brightest source in this sub-cluster 
is labeled IRS3. 



The fourth cluster is associated with IRAS 16177-5018 
and, together with IRAS 16132-5039, is embedded in the Hn 
region RCW 106. Roman- Lopes et al. (2003), obtained JHK 
images and photometry of the reddest stars, which have 
visual extinction exceeding 28 magnitudes. The brightest 
source is IRS1, located at the centre of the infrared nebula, 
which together with IRS7, seem to play a key role in the en- 
ergy balance of the compact Hn region (Roman-Lopes et al. 
2003). An image of the region, with the two stars indicated, 
is shown in Figure Id. 



3 OBSERVATIONS AND RESULTS 
3.1 GNIRS data 

7f-band spectra of the point sources indicated in Figure 1 
were obtained on different dates using GNIRS (Elias et al. 
2006) on the 8-m Gemini South telescope at Cerro Pachon, 
Chile. Table 1 shows the log of observations. In all cases the 
short camera with a pixel scale of 0.15" pixel -1 and a slit 
size 0.675" x 99" was employed. The resolving power of this 
configuration is 1600, with a theoretical wavelength coverage 
ranging from 1.90 to 2.50 /im. 

In order to obtain good results in the sky subtraction 
process, the standard ABBA nodding technique was used to 
acquire the spectra. The large slit length (99") enabled us 
to perform large node shifts with the target still on the slit. 
This is specially useful when observing regions where the 
nebular lines (like Br7) are spatially extended, to avoid the 
overlap of AB positions that would create artifacts in the 
sky subtracted images. Individual exposure times at each 
nod position was 3 minutes for the science targets; total 
exposure times are listed in column 5 of Table 1. In order 
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Figure 4. GNIRS K-band spectra of the hot stars detected in this work. All spectra were flux normalized and are at the same scale. 



to correct the science spectra for the effect of telluric at- 
mospheric absorption, a nearby AOv spectroscopic standard 
star was observed at similar airmass (column 7 of Table 1), 
before or after the set of exposures of each science target. 

The GNIRS data were reduced using the GEMINI pack- 
age within iraeQ. First, the 2D A"-band frames were sky- 
subtracted for each pair of images taken at the two nod po- 
sitions A and B, followed by division of the resultant image 
by a flat-field. Multiple exposures for each source were com- 
bined, followed by one-dimensional extraction of the spectra. 



IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, INC, under cooperative agreement with 
the National Science Foundation 



Eventually, wavelength calibration was performed using sky 
lines; the typical error (1-cr) for this calibration was ~5 A. 
Telluric atmospheric correction using the spectroscopic stan- 
dard stars completed the reduction process. In this last step, 
we divided the target spectra by the spectrum of the AOv 
spectroscopic standard star, already free of photospheric fea- 
tures. In the standard star, the Br7 absorption is the only 
feature present in the if-band spectrum. It was carefully 
removed by interpolation across its wings using continuum 
points on either side of the line while its core was modeled 
using a Voigt profile. In order to assure good cancellation 
of the telluric bands, the IRAF task telluric was employed. 
The algorithm interactively minimizes the RMS in specified 
wavelength regions by shifting and scaling the target rel- 
ative to the standard spectrum to best divide out telluric 
features present in the former. Shifting account for small er- 
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rors in the dispersion zero-points, while the intensity scaling 
corrects for differences of airmass and variations in the abun- 
dance of the telluric species. Typical values of the shifts were 
a few tenths of a pixel (equivalent to ~2 A) and the scaling 
factors were less than 10%. As an example, Figure 2 shows 
the spectra of the AOv standard star HIP83818 (already free 
of the Br7 line), the science target IRAS16177-5018JRS7, 
and its spectrum corrected from the effects of the telluric 
absorption bands. 

3.2 Correcting the science spectra for nebular 
contamination 

As can be inferred from Figures lb and Id, the re- 
gions associated to IRAS 15408-5356 and IRAS16177-5018 
present strong extended emission. In fact, the K-h&nd spec- 
tra of IRAS15408-5356JRS1, IRAS16177-5018.IRS1, and 
IRAS16177-5018JRS7 are contaminated by Br7 and Hel 
(2.058 pirn) nebular components. Since these lines play a 
fundamental role in the spectral classification, they must be 
carefully subtracted from the stellar spectrum. To do that, 
we evaluated the nebular contribution at the position of the 
point sources by studying the intensity profiles of the Br7 
and Hel (2.058 fim) lines as function of the position along 
the slit. 

In Figure 3 we illustrate the procedure used in the case 
of the sources belonging to the IRAS16177-5018 cluster. At 
the bottom panel of this figure, we present the contour di- 
agram made from the GNIRS acquisition image; there we 
indicate the position on the slit of IRS1 and IRS7, with 
the spatial coordinate scale measured in pixels. Figures 3a 
and 3d show the spectral intensity at the wavelength of the 
Br7 and Hel 2.058 /im lines respectively, which include the 
continuum emission. After the subtraction of the continuum, 
taken at ±30 A(about 6 pixels) off the line centre and shown 
in Figures 3b and 3e, we obtained the Br7 and Hel 2.058 
/im line contributions, presented in Figures 3c and 3f, re- 
spectively, where we can see the nebular emission at both 
sides of the stars. In Figure 3f Hel 2.058 /jm is clearly seen 
in absorption at the position of IRS1, while no absorption 
is present in Br7 in Figure 3c. This procedure was used to 
obtain the spectra of all stars embedded in ionised clouds, 
even when the nebular contribution was small. 



3.3 Spectral classification 

The GNIRS spectra have been organized in three groups. 
In the first one, the detected lines are characteristic of hot 
stars (Br7, Hel 2.058 fim, Hel 2.113 //m, Hell 2.185 //m, Civ 
2.078 /jm, Niii 2.116 fim, among others.); they are shown in 
Figure 4 in order of increasing Br7 equivalent width (EW). 

The spectra of the sources in the second group, shown 
in Figure 5, present weak Br7, metallic (Na, Ca, Mg, etc.), 
and molecular (CO overtones) lines, characteristic of late- 
type stellar spectra. Finally the third group shows H and 
He lines in emission but no photospheric absorption lines, 
which is characteristic of young stellar objects (YSOs); their 
spectra are shown in Figure 6. 

The library of K-bwd early-type spectra compiled by 
Hanson, Conti & Rieke (1996) was used to classify the stars 
in our sample. Since both spectra have similar resolution, 
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Figure 5. GNIRS spectra of the late-type stars detected in this 
work. All spectra were continuum flux normalized and are at same 
scale. 



the classification of the stars in group 1 was made by direct 
comparison. For each star, we show the library spectrum 
that best matches the target one, as well as earlier and later 
type library spectra (Figures 7, 8, 9 and 11). 

The derived spectral types are presented in Table 2, 
together with the absolute Mj magnitudes. They were com- 
puted considering the Mv magnitudes taken from Hanson 
ct al. (1997) and Walborn (2002), and the intrinsic colours 
given by Koornneef (1983) transformed into the 2MASS 
photometric system. Columns 6 and 7 list the intrinsic 
(J — H)o and (H — K )o colours corresponding to each spec- 
tral type. Columns 8, 9, 10 and 11 list the measured J and 
K magnitudes, as well as the (J — H) and (H — K) colours 
obtained from previous works (Roman-Lopes et al. 2003, 
Roman-Lopes & Abraham 2004a, 2004b, Ortiz et al. 2007). 



3.4 Determination of the distances 

The distances d to the stars were determined using abso- 
lute and measured magnitudes presented in Table 2, in the 
equation: 

rn x — M\ = 5 log[d(pc)] - 5 + A(\) (1) 

The absorption A(X) is related to the colour excesses 
E(J - H) and E(H - K) through the functions Fj(R) = 
A(J)/E(J - H) and F K (R) = A(K)/E(H - K) that de- 
pend on the ratio of the total to selective extinction R — 
Av /E(B — V). Using the interstellar extinction laws given 
by Fitzpatrick (1999), it is possible to obtain the ratios 
Aj/A v = fj(R), Ak/Av = f K (R), E(J - H)/E(B -V) = 
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Figure 6. GNIRS spectra of YSO's detected in this work. All 
spectra were continuum flux normalized. 



gj(R) and E(H - K)/E(B - V) = g K {R) from which the 
functions F(R) = Rf(R)/g(R) can be derived. 

A source of uncertainty in the distance determination 
originates from the use of the standard interstellar extinction 
law, which is represented by R — 3.1, though it is common 
to find R in the range 2.8 — 5.8 along the galactic plane 
(Johnson 1965, Tapia 1981, Fitzpatrick 1999, Indebetouw 
et al. 2005, Nishiyama et al. 2006). This effect is probably 
produced by differences in metallicity and grain size distri- 
bution (Savage & Mathis 1979). It was taken into account 
by using R = 2.8, 3.1 and 5.0, and the distance to each star 
was calculated as the average of these individual values. 

Another important source of error comes from the spec- 
tral type determination itself, which in our work has an un- 
certainty of about ± one sub-type. For O-type stars it repre- 
sents about 0.2 magnitudes in the NIR, whereas for early-B 
stars it is about 0.6 magnitude (Hanson et al. 1997). The 



only exception is the I16177-5018_IRS1 source, for which the 
true luminosity class is an important source of uncertainty, 
this issue will be discussed in detail in Section 4. 

A third source of uncertainty results from the fact that 
young stars can be surrounded by disk and/or dust cocoons, 
which produce excess emission in the NIR (Grasdalen et al. 
1975, Glass 1979, Lada & Adams 1992), especially in the 
K-band. In order to minimize this effect, we used both the 
J and K band photometry to derive distances. 

Table 3 shows the values of E(J - H), E(H - K), A,, 
and Ak for three values of total to selective extinction ratios 
(2.8, 3.1 and 5.0). The stellar distances d were computed as 
the average of the distances obtained from the J and K mag- 
nitudes, and the three extinction laws. The quoted errors for 
the distances are the standard deviation of the averaged val- 
ues, which include the errors associated to the photometry, 
interstellar extinction law, and spectral type, as discussed 
above. Finally, the distance to each cluster (Column 13) was 
computed as the average of the distances to each star in the 
cluster; the total error of these values is the square root of 
the sum in quadrature of the individual errors. 



In this section, we discuss our main results and compare 
them with those obtained in the previous works. 



4.1 The IRAS09149-4743 cluster 

In this region, three program stars were chosen for observa- 
tion: IRS1, IRS2 and IRS3. In addition, three other sources 
located nearby had their spectra taken as well: IRS18, 
IRS25, and IRS31 (Ortiz et al. 2007). Their relative posi- 
tions are shown in Fig. la. The source labeled as IRS1 was 
resolved into two stars in this work: IRS la and IRS lb (see 
the inset in Fig. la). The spectra of IRSla and IRS2, shown 
in Fig. 4, exhibit two important diagnostic lines commonly 
found in stars earlier than Blv: Hi Br-7 at 2.166 /im and 
the Hel line at 2.113 /im. The source IRSla was classified as 
a BOv star (Figure 7a), affected by 9 magnitudes of visual 
extinction (Table 3), which places it at a distance of about 
1.20 ±0.12 kpc. 

IRS 2, belonging to the nearby sub-cluster located 1.1 
arcmin to the southeast of IRS1, is probably an 09v star 
(Figure 7b) affected by about 7—8 visual magnitudes; it is 
located at 1.27 ± 0.13 kpc, virtually the same distance as 
IRSla. IRS31 belongs to the same sub-cluster as IRS2. Its 
spectrum (Fig. 7c) exhibits a strong Br-7 line, characteristic 
of late-B and early-A stars. This source has been classified 
as a B7v-B8v star, suffering 5-8 magnitudes of extinction in 
the V band (Table 3). Its distance is about 1.37 ± 0.19 kpc, 
reinforcing its membership. 

The distances to the three massive stars studied in this 
region are similar, resulting in a mean cluster distance of 
1.3 ± 0.2 kpc. This value can be compared with the pho- 
tometric distance of 0.7 ± 0.2 kpc quoted by Liseau et al. 
(1992) for cloud A in the VMR complex, and the kinematic 
distance of 1 kpc inferred from CO observations by Murphy 
& May (1991). The good agreement between our result and 
that obtained from the Galactic rotation curve is notable, 
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Figure 7. GNIRS K band spectra of IRAS09149-4743 sources 
(red lines) overlayed on iv-band spectra (black lines) taken from 
the library of Hanson et al. (1996). 



considering the complexity of the VMR and the small radial 
velocity resulting from its position close to I = 270° . 

In Figure 6 we can see that IRS lb shows no photo- 
spheric spectral lines, but Br7 in emission and some CO 
overtone band-heads in absorption, characteristic of YSO's 
(Casali & Eiroa 1996, Hoffmeister et al. 2006). These CO 
lines are believed to be formed in a warm and dense circum- 
stellar shell, possibly a relic of a former accretion disk. They 
can be seen in absorption or emission depending on the disc 
opacity, which in turn depends on the mass accretion rate, 
as shown by Calvet et al. (1991). Ortiz et al. (2007) also 
pointed out that IRSla+b shows intense infrared emission 
beyond 5 /im, usually attributed to warm dust. Based on the 
present data, we can now state that the observed infrared 
excess comes from IRS lb, the dusty nearby companion of 
the IRSla source. 

IRS3 is a highly reddened IR source (Ortiz et al. 2007), 
located at about 1 arcmin west of IRS1 (Figure la). Its spec- 
trum, shown in Figure 5, exhibits metallic and molecular 
absorption lines, such as the CO (2,0) and (3,1) transitions 
at 2.29 and 2.32 /im respectively, as well as the Cal and Nal 
lines at 2.21 and 2.26 /im, typical of K and M giant stars. As- 
suming its spectral classification as early-K, we can estimate 
a lower limit for its distance. Its absolute magnitude and in- 
trinsic colour index would be Mj = —1.9 and (J—H) = 0.62, 
which implies a colour excess E(J — H) — 1.42. If one as- 
sumes R = 3.1 then Aj — 4.34 and d = 2.5 kpc. On the 
other hand, if R — 5.0, Aj — 4.76, and d — 2.0 kpc. In any 
case, the distance to this star would be twice as large as the 
distance to RCW 41. Besides IRS3, two other stars in the 
neighbourhood have been classified as late- type: IRS 18 and 
IRS25 (Fig. 5). Similarly to IRS3, they occupy a position in 
the (J — H) versus (H — K) C-C diagram consistent with 
their classification as late-type stars, and their JHK magni- 
tudes given in Ortiz et al. (2007) imply that they must also 
be background objects. 

4.2 The IRAS15408-5356 cluster 

In this cluster, spectra of four program stars (IRS1, IRS2, 
IRS3, and IRS4) were obtained. The spectrum of IRS1 is 
shown in Figure 4. It presents features typical of very hot 
stars, such as Civ (at 2.078 /jm) and Nin (at 2.116 fim) 
in emission, characteristic of stars with spectral type earlier 
than 06v. A direct comparison with the library templates 
of Hanson et al. (1996) allowed us to classify it as 05. 5v 
(Figure 8a). Morisset (2004) found that the relative intensi- 
ties of mid-infrared emission lines observed by the Infrared 
Space Observatory (ISO), require an ionizing source with ef- 
fective temperature T e ff = 48, 700 K, equivalent to an 03v 
star, implying that other stars must be contributing to the 
total ionizing luminosity. 

IRS3, similar to IRS1, was previously suggested to be 
one of the main ionising sources of RCW 95 (Roman-Lopes 
& Abraham 2004b). It actually consists of two stars, labeled 
IRS3a and IRS3b, as shown in the inset in Figure lb. The 
former presents an unexpected spectrum; it does show the 
hydrogen Br7 and the Hel lines typical of hot stars, but also 
the CO band-head overtone lines in absorption, characteris- 
tic of late-type stellar atmospheres. Although CO overtone 
bands in absorption have been widely reported in low-mass 
YSO spectra (Straw, Hyland & Mc Gregor 1987, Straw et 
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Figure 8. GNIRS spectra of IRAS15408-5356 IRS1, IRS3a and 
IRS3b sources (red lines) overlayed on K band spectra (black 
lines) taken from the library of Hanson et al. (1996). 



al. 1987, Carr 1989, Casali & Matthews 1992), they have 
also been found in a few late-O/early-B stars (Hoffmeister 
et al. 2006). The origin of these CO features is not clear; 
they might be the signature of a cold star in the same line 
of sight or simply neutral gas in the interface between the 
HII region and the molecular cloud. Apart from the CO lines 
beyond 2.3 /im, this star can be classified as 09. 5v, as shown 
in Figure 8b. IRS3b has been classified as an B3/B5v star 
(Figure 8c). 

Bik et al. (2005) obtained isT-band spectra for two stars 
in common with the studied region: 15408nrl410 (05v- 
06.5v) and 15408nrl454 (08v-B2.5v) and, according to 
their coordinates (Bik 2004), these objects correspond to 
sources IRS1 (05.5v) and IRS3a (09. 5v). Therefore, the 
spectral classification of the two works agree with each other 
within the error bars. 

The spectrum of IRS2, seen in Figure 6, does not 
present any evident photospheric features, except for the 
Br7 emission line, superposed on an absorption profile. An- 
other emission line near 2.27 /im might be due to the 4- 
3 S(4) H 2 transition. Differently from IRAS 09149-4743 
IRSlb, the spectrum of IRS2 does not show any evidence 
of a CO environment that could be associated with circum- 
stellar disks. However, since this object shows large IR color 
excess and is associated with the bright source MSX6C- 
G326. 6570+00. 5912, our result confirms its previous clas- 
sification as an YSO (Roman-Lopes & Abraham 2004b). 

Three additional less luminous stars located near IRS 2 
fell into the slit: IRS 10 and other two sources not included 
in the previous photometric study of the region, labeled 
IRS2_ap2 and IRS2_ap3. Their spectra, shown in Figure 5, 
are characteristic of late-type stars. This result implies that 
IRS10 is not one of the ionising sources of RCW95, as pre- 
viously proposed considering its JHK magnitudes and col- 
ors (Roman-Lopes & Abraham 2004b). Its location in the 
colour-colour diagram seems to result from its low effective 
temperature and the high extinction of the cloud, which to- 
gether mimic the colour indices of young, massive stars. 

The spectrum of IRS4, shown in Figure 5, is clearly 
late-type, with strong CO absorption lines in the spectral 
region between 2.29 /im < A < 2.4/xm. One can see also a 
few metallic lines, probably due to Nal and Cal transitions, 
which reinforces this classification. 

The distances to the early-type sources have been de- 
termined from their spectral types, J and K magnitudes, 
and E(J — H) and E(H — K) colour excesses, as described 
in Section 3.4, and presented in Table 3, with the excep- 
tion of source IRS3b. Since no photometric data are avail- 
able for this object, we used the visual absorption calcu- 
lated for IRS3a instead, and the K magnitude obtained 
from the GNIRS A"-band acquisition image. The instrumen- 
tal GNIRS magnitudes were calibrated using the combined 
K magnitudes of IRS3a and IRS3b, obtained from the work 
of Roman-Lopes & Abraham (2004b). The resulting indi- 
vidual K magnitudes are 9.5 and 12.0 for IRS3a and IRS3b, 
respectively. 

The derived distances to all stars observed in this cluster 
are similar, like in IRAS09149-4743, ranging from 1.32 ±0.20 
kpc for IRS3a to 1.34±0.16 kpc for IRS1, giving a mean dis- 
tance of 1.3±0.2 kpc for the IRAS 15408-5356 cluster. Giveon 
et al. (2002) derived a kinematic distance of 2.4 kpc using the 
Galactic rotation curve, assuming galactocentric distance of 
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8.5 kpc and a solar rotation velocity of 220 km s . This 
is about twice what we found, showing that at least in this 
direction, the Galactic kinematic model fails. In fact, differ- 
ences between the kinematic distances, derived from radio 
observations and those obtained by spectroscopic paralaxes 
were also found by Blum et al. (1999, 2000) and Figueredo 
et al. (2005), showing the importance of spectrophotometric 
studies for a better understanding of the rotation curve of 
the Galaxy. 

4.3 The clusters in the RCW106 region 

The two clusters associated to IRAS16132-5039 and 
IRAS16177-5018 are part of the RCW106 complex. One 
program target was observed in each star formation region: 
IRS3 in the former cluster and IRS1 in the latter. 

The spectrum of the IRS1 source (Fig. 4), shows Civ, 
Nni, and Nv emission lines at 2.078 /jm, 2.116 /jm, and 
2.100 fim, respectively, Hell at 2.189 /im in absorption, and 
a weak Br7 line. The presence of the Nv line indicate that 
IRS1 might be a very hot star. In fact, such line is only found 
in the A'-band spectra of 03v and 03-04lf + supergiant stars 
(Hanson et al. 2005). 

In Figure 9a we compare the IRS1 spectrum with that 
obtained for OS-OSlf 1 " supergiant stars (Hanson et al. 1996). 
One can notice that the GNIRS spectrum resembles well 
that of Cyg OB2 7, an OSlf 1 " star. For completeness, in Fig. 
9b we also compare the IRS1 spectrum with that of 03- 
06 main-sequence stars. In this case we see a reasonable 
agreement between the GNIRS spectrum features with that 
from the HD93250 (03v), HD164794 (04v), and HD93204 
(05v). The exception are the Br7 and He n lines, which in 
the templates appear stronger in absorption. 

There are additional constraints indicating that the 
main ionizing source of the compact Hn region probably 
is an extremely hot source. From the K-band spectrum of 
the associated nebular emission, which is shown in Fig. 10, 
we can see the presence of strong emission lines like the 
Br 7 (2.166 fim), Hel (2.058 fim, 2.113 fim and 2.161-2.162 
pirn), and [Fein], which are normally found in HII regions 
hosting embedded hot stars. We also found the s-process 
[Krin] and [Seiv] high excitation emission lines, previously 
identified only in planetary nebula (Sterling et al. 2007). 
Blum & McGregor (2008) also detected such lines in their 
study of the ionising stars associated to the UC Hn region 
G45.45+0.06, suggesting that high density Hn regions ex- 
cited by the hottest O stars also would produce such emis- 
sion lines. Indeed, the measured Hei 2.113/Br7 ratio (see 
Fig. 10) of 0.045±0.003 (not corrected for reddening) indi- 
cate the presence of an exciting star whose T e fj is greater 
than 40,000 K (Hanson, Luhman & Rieke 2002). This would 
correspond to a star of spectral type earlier than 05v, 
04.5III, and 04l (taking into account the luminosity class), 
as inferred from the new calibration of O star parameters 
published by Martins et al. (2005). 

The fact that IRS1 has a K-band spectrum similar to 
that of Cyg OB2 7, the detection of the s-process [Krin] and 
[Seiv] high excitation emission lines (from the associated 
UC Hn region), and the lower limit to the effective temper- 
ature (40,000 K) for the main ionising source in the region, 
is reasonable to suppose that IRS1 may be a 03 supergiant 
star (or alternatively a very hot ionising source). Taking into 
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account the rarity of such type of objects, this is an extraor- 
dinary result, which could indicate that RCW106 may be 
the birthplace of extremely massive stars. 

The presence of a supergiant star in a young massive 
stellar cluster could be questioned in view of the youth of 
the sources in the region (Roman-Lopes, Abraham & Lepine 
2003). In fact, Roman- Lopes (2007) using the fraction of 
NIR sources showing excess emission in the NIR, estimated 
the age of the cluster as 2.5-3.0 Myrs. Indeed, other O-type 
giants and supergiants have been found in very young stellar 
clusters. For instance, Massey et al. (2001) derived ages for 
several massive early-type stars, founding values in the range 
1.0-3.0 x 10 6 years. Recently, Melena et al. (2008) detected 
several early O giant and supergiant stars with similar ages. 
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Figure 10. (a) GNIRS K band spectra of the 116177 UC HII 
region. The intensity of the lines were normalized to that of the 
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excitation emission lines. 



Assuming that IRS1 is a 03 If+ supergiant star, we 
computed its distance as 2.6±0.7 kpc, compatible with that 
of 2.8±0.6 kpc obtained from the galactic rotation curve 
assuming Ro = 8.5 kpc, 6 = 220 km s" 1 (Honma & Sofuc 
1997), and a radial velocity Vr=-49.5 km s _1 . On the other 
hand, considering IRS1 as a 03-05 v star, the distance drops 
to about 1.2±0.7 kpc. Indeed, further 7f-band spectroscopic 
observations of other nearby NIR sources will be necessary 
to improve our understanding of the Galactic structure in 
this direction. 

Bik et al. (2005) also used NIR spectroscopy to estimate 
the distance to this cluster. Three sources were observed 
(16177nr271, 16177nr405, and 16177nrl020), but unfortu- 
nately there are no coordinates or other information that 
would allow us to perform a cross-correlation between those 
data with ours. The individual determinations of the dis- 
tance in that work show somewhat large error bars (1.0 — 
3.7 kpc), but considering its mean value of about 2.4 kpc, 
their results can be considered consistent with ours. 

Another observed point source in this cluster is IRS7 
which spectrum does not show any photospheric spectral 
features but the Br7 and Hel (2.059 /im) lines in emission, 
characteristics of an YSO (Figure 6). As pointed out by 
Roman-Lopes et al. (2003) this object coincides with the 
mid-infrared source MSX6C-G333.3072-00.3666, which re- 
inforces this classification. We notice a P-Cygni profile in 
the Hel line at 2.058 /im that may indicate the presence of 



o 
+ 



B 

u 

o 



Hel Br-y 
l He! ' 



HD 38666 



09.5 V 



HD 36512 



B0 V 



HD 36960 B0.5 V 



! 



_L 



2,04 3.08 2.12 2.16 
Wavelength (fJ,m) 

IRAS16132-5039_IRS3_ap3 



2.20 



u 

Xfi 

o 



tn 

01 



t3 
S 

o 
2 



1 

Hel 


Br-y 

Hel I 
! 


_HD|£1071 


B7 V . 


HDj\169990 


V B8 V 


HD .97583 


B8/9 V 


T 

1 

1 1 1 



2.04 2.08 2.12 2.16 
Wavelength (^m) 



2.20 



Figure 11. GNIRS K band spectra (red lines) of two sources in 
the cluster associated to the IRAS16132-5039 source, overlayed 
on K band spectra (black lines) taken from the library of Hanson 
et al. (1996). 



expanding material inside the region. High-resolution radio 
observations would be useful to clarify this issue. 

The other studied cluster in RCW106 is IRAS16132- 
5039. During the observation of IRAS16132-5039JRS3 three 
stars fell into the slit: IRS3 itself, which has a K band spec- 
trum similar to O9.5v/B0v (Figure 10a), and two additional 
sources: IRS3_ap2, identified as an YSO (Figure 5), and 
IRS3_ap3, probably a B8V/B9V star (Figure 10b). Similarly 
to IRAS15408-5356JRS3a, the spectrum of IRAS16132- 
5039JRS3 presents H and He lines as well as CO overtone 
lines; the possible origin of the CO lines was already dis- 
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cussed in a previous section. The distances to IRS3 and 
IRS3_ap3 were found to be 2.15 ± 0.50 and 2.64 ± 0.30 kpc, 
respectively (Table 3). Therefore, it can be concluded that 
both sources are located at the same distance. The derived 
spectrophotometric mean distance is about 2.4±0.5 kpc, a 
value that is compatible with that obtained from the Galac- 
tic rotation curve (3.2±0.6 kpc) using the radial velocities 
taken from the CS (2-1) line (Bronfmann et al. 1996) and 
from the hydrogen recombination lines (Caswell & Haynes, 
1987). 



5 SUMMARY 

We report A"-band spectroscopic observations of stars 
in highly reddened young stellar clusters, obtained with 
GNIRS. We found 8 massive main-sequence stars, 1 possible 
supergiant, 7 late-type stars and 4 YSOs. 

The main-ionising star in the cluster associated to 
IRAS09149-4743, is IRS2 (09v), which together with IRS1 
(BOv) are the dominant ionising sources in the RCW41 HII 
region. They are located in cloud A of the VMR at a mean 
distance of 1.3±0.2 kpc, in agreement with the value derived 
from the kinematic method using CO radial velocities. We 
found that IRS09149JRS1 actually consists of two sources: a 
BOv star (IRSla) and an YSO (IRSlb). This result allows us 
to state that the infrared excess found by Ortiz et al. (2007) , 
comes from the dusty nearby companion of IRSla. The other 
bright source in the region, IRS09149 JRS3, is a background 
giant star located at least twice as far as the cluster, con- 
firming the assumption made by Ortiz et al. (2007), who 
suggested that this source is a late-type star. 

Four A'-band spectra were obtained in the cluster asso- 
ciated with RCW95. IRAS15408JRS1 has a if-band spec- 
trum compatible with an 05. 5v star, while IRS3a and IRS3b 
were classified as 09. 5v and B3-B5v, respectively. Their dis- 
tances are similar and have an average value of 1.3±0.2 kpc, 
in disagreement with the kinematic distance of 2.4 kpc. This 
seems to indicate that, at least for this star forming region, 
the galactic kinematic model fails. Eventually, these results 
can be used to improve our understanding of the rotation 
curve of the Galaxy. IRS2, another very bright source in 
this cluster, is associated with the mid-IR source MSX6C- 
G326. 6570+0. 5912, and was classified by Roman-Lopes & 
Abraham (2004b) as an YSO candidate. Its if-band spec- 
trum presents Br7 and Hel emission lines characteristic of 
YSOs, confirming the previous classification. 

One of the sources observed in the direction of 
RCW106 belongs to the cluster associated with IRAS16132- 
5039. Three other stars fell into the spectrograph slit: 
IRAS16132JRS3 (O9.5-B0v), IRAS16132_IRS3.ap2 (YSO) 
and IRAS16132JRS3_ap3 (B8-B9v). The derived mean dis- 
tance is 2.4±0.5 kpc, a value compatible with that obtained 
from the kinematic method (3.2±0.6 kpc) using the CS (2-1) 
and hydrogen recombination lines. 

The other cluster in the RCW106 region is associ- 
ated with IRAS16177-5018. The IRAS16177JRS1 spectrum 
presents Civ, Niii, and Nv emission lines at 2.078 /jm, 2.116 
/jm, and 2.100 /im respectively, which are seen only in very 
hot stars. We also detected in the associated nebular K- 
band spectrum, the s-process [Krlll] and [Seiv] high excita- 
tion emission lines, previously identified only in planetary 



nebula. In the case of Hn regions, such lines seem to be 
produced only in high density environment excited by the 
hottest O stars. Indeed, the measured Hel 2.113/Br7 ra- 
tio indicate the presence of an exciting star whose T e // is 
greater than 40,000 K. This would correspond to a single 
star of spectral type earlier than 05v, or 04.5ni, or 04l, 
depending on the luminosity class. 

We found a good agreement between the 
IRAS16177JRS1 if-band spectrum and that of Cyg 
OB2 7, an 03If* star. Taking into account the scarcity 
of supergiants in young cluster (though a few similar 
occurrences have been reported elsewhere), this is an 
extraordinary result, which could indicate that RCW106 
may be the birthplace of extremely massive stars. 

Considering IRAS16177JRS1 as an 03If* star, the dis- 
tance to the cluster would be 2.6±0.7 kpc, similar to that 
of the other massive star formation region in the RCW106 
complex. On the other hand, if IRAS16177JRS1 is an 03- 
05 main-sequence star, its distance drops to 1.2±0.7 kpc, 
much smaller than its kinematic distance. 

We originally planned to obtain the K-band spectrum 
for one target in the IRAS16177-5018 region, but another 
source also fell into the slit: IRAS16177JLRS7. Its spectrum 
is typical of an YSO, in agreement with the classification 
proposed by Roman-Lopes, Abraham & Lepine (2003). Its 
Hel (2.058 /im) emission line shows a P-Cygni profile that 
may indicate the presence of expanding motion. High resolu- 
tion radio observations would be useful to clarify this issue. 
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Table 1. Log-book of the GNIRS observations. The columns are: (1) The associated IRAS source; (2) the identifier of the 
point source, as designated in the original photometric survey;(3,4) equatorial coordinates (J2000.0);(5) the total exposure 
time (s);(6)the mean airmass at the time of the observations; (7) the Hipparcos identifier of the associated AOV telluric 
star; (8) the mean airmass of the telluric star at the time of the observations; (9) lists the Gemini identification program and 
Column (10) the corresponding date of observation. 



Source 


IRS 


a (J2000.0) 8 (J2000.0) 


Itime (s) 


X Telluric 


X 


ID Program 


Date 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) (7) 


(8) 


(9) 




(10) 




09149-4743 


1 


09hl6m43.50s -47°56'23.0' 


' 540 


1.07 HIP40974 


1.03 


GS-2005B- 


Q-33 


Dec 23, 


2005 




09149-4743 


2 


09hl6m47.94s -47°57'18.0' 


132 


1.06 HIP40974 


1.05 


GS-2005B- 


Q-33 


Dec 23, 


2005 




09149-4743 


3 


09hl6m41.89s -47°56'15.4' 


540 


1.05 HIP40974 


1.06 


GS-2005B- 


-Q-33 


Dec 24, 2005 




15408-5356 


1 


15h44m43.40s -54°05'53.7' 


' 180 


1.30 HIP75161 


1.36 


GS-2005B- 


Q-33 


Dec 24, 2005 




15408-5356 


2 


15h44m43.40s -54°05'53.7' 


180 


1.34 HIP75161 


1.35 


GS-2005B- 


Q-33 


Dec 27, 2005 




15408-5356 


3 


15h44m56.17s -54°07'18.1' 


360 


1.57 HIP75161 


1.63 


GS-2005B- 


■Q-33 


Jan 22, 2006 




15408-5356 


4 


15h44m53.04s -54°06'31.9' 


' 1080 


1.64 HIP75161 


1.63 


GS-2005B- 


■Q-33 


Jan 23, 2006 




16132-5039 


3 


16hl6m55.99s -50°47'22.8' 


1080 


1.88 HIP83406 


1.75 


GS-2005B- 


Q-43 


Sep 25, 2004 




16177-5018 


1 


16h21m31.60s -50°25'08.3' 


1080 


1.56 HIP83818 


1.65 


GS-2005B- 


Q-33 


Jan 24, 


2006 




Table 2. Derived spectral types for the hot stars in our 


sample. 
















Source 


a 


(J2000) 


8 (J2000) 


ST. 


Mj 


(J-H) Q 


(H-K) 


J 


(J-H) 


K 


(H-K) 


(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


I09149JRSla 


09hl6m43.50s 


-47°56'23.0" 


B0V 


-2.21 


-0.16 


-0.01 


10.70 


0.80 


9.40 


0.50 


I09149.IRS2 


09hl6m47.94s 


-47°57'18.0" 


09V 


-3.00 


-0.18 


-0.01 


9.72 


0.66 


8.68 


0.38 


I09149JRS31 


09hl6m47.77s 


-47°57'15.1" 


B7/B8V 


-0.27/0.29 


-0.07 


+0.02 


13.34 


0.80 


12.25 


0.29 


I15408.IRS1 


15h44m43.40s 


-54°05'53.7" 


05.5V 


-4.05 


-0.20 


-0.01 


10.76 


1.37 


8.57 


0.82 


I15408JRS3a 


15h44m56.17s 


-54°07'18.1" 


09.5V 


-2.85 


-0.16 


-0.01 


11.25 


1.17 


9.50 


0.58 


I15408.IRS3b 1 


15h44m56.17s 


-54°07'18.1" 


B3/B5V 


-0.45/+0.15 


-0.11 


+0.01 






12.0 




I16132_IRS3 


16hl6m56.01s 


-50°47'22.7" 


O9.5/B0V 


-2.85/-2.21 -0.18/-0.16 


-0.01 


11.77 


0.76 


10.30 


0.71 


I16132_IRS3.ap3 


16hl6m50.84s 


_50°47'44.2" 


B8/B9V 


-0.27/0.29 


-0.03 


-0.01 


12.71 


0.15 


12.50 


0.06 


I16177.IRS1 


16h21m31.60s 


-50°25'08.3" Q3If*/03-05V 


-6.85/-4.95 


-0.21 


-0.01 


16.65 


3.92 


10.21 


2.52 



Table 3. Summary of the extinction values and distances to the program stars. For the I16177JR.S1 source we show the distances obtained 
considering class I (1) and class V (2) case. 



Source 


E(J - H) 


Aj 


Aj 


Aj 


A v 


E(H - K) 


A K 


A K 


A K 


A v 


d 


clusi 






R = 2.8 


R = 3.1 


R = 5.0 


R = 3.1 




R = 2.8 


R = 3.1 


R = 5.0 


R = 3.1 


(kpc) 


(kpc 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13 


I09149JRSla 


0.96 


2.5 


2.6 


2.8 


9.1 


0.51 


1.0 


1.0 


1.1 


8.9 


1.20 + 0.12 


1.3 + 


I09149.IRS2 


0.84 


2.2 


2.3 


2.5 


8.0 


0.39 


0.9 


0.9 


1.0 


6.8 


1.27 + 0.13 




I09149.IRS31 


0.87 


2.3 


2.3 


2.6 


8.3 


0.27 


0.9 


0.9 


1.0 


4.7 


1.37 + 0.19 




I15408JRS1 


1.57 


4.1 


4.2 


4.6 


14.9 


0.83 


1.6 


1.7 


1.9 


14.4 


1.34 + 0.16 


1.3 + 


I15408JRS3a 


1.33 


3.5 


3.6 


3.9 


12.6 


0.59 


1.4 


1.4 


1.6 


10.2 


1.32 ± 0.20 




I15408JRS3b 


* 


* 


* 


* 


12.6 


* 




* 


* 


10.2 


1.24 + 0.21 




I16132JRS3 


0.92 


2.4 


2.5 


2.7 


8.7 


0.72 


1.0 


1.0 


1.1 


12.5 


2.15 ± 0.50 


2.4 + 


I16132.3(ap3) 


0.22 


0.6 


0.6 


0.7 


2.1 


0.04 


0.2 


0.2 


0.3 


0.7 


2.64 ± 0.30 




116177 JRSlW 


4.13 


10.7 


11.1 


12.1 


39.2 


2.53 


4.3 


4.5 


4.9 


43.9 


2.58 ± 0.68 


2.6 ± 


I16177JRS1< 2 > 


4.13 


10.7 


11.1 


12.1 


39.2 


2.53 


4.3 


4.5 


4.9 


43.9 


1.22 ± 0.74 


1.2 + 



